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We discuss finite-size effects on homogeneous nucleation 
in first-order phase transitions. We study their implications 
for cosmological phase transitions and to the hadronization 
of a quark-gluon plasma generated in high-energy heavy ion 
collisions. Very general arguments allow us to show that the 
finite size of the early universe has virtually no relevance in the 
process of nucleation and in the growth of cosmological bub- 
bles during the primordial quark-hadron and the electroweak 
phase transitions. In the case of high-energy heavy ion colli- 
sions, finite-size effects play an inportant role in the late-stage 
growth of hadronic bubbles. 



I. INTRODUCTION 

First-order phase transitions and the kinetic phenom- 
ena associated with the process of phase conversion 
through the nucleation of bubbles or spinodal decom- 
position are present in almost all realms of physics [1]. 
Usually, the interesting physical quantities are the time 
scales related to the transient regime, since they deter- 
mine its relative importance. However, one cannot often 
disentangle those scales from other scales of the problem, 
such as the different length scales and, in particular, the 
finite size of the system under consideration. 

Despite the success of finite-size scaling in the study of 
equilibrium critical phenomena [2] , a systematic study of 
finite size effects has rarely been conducted in the case 
of metastable decays and other nonequilibrium processes 
(see Ref. [3] and references therein). In particular, a care- 
ful analysis of possible modifications of the usual picture 
of nucleation in a thermally-driven first-order phase tran- 
sition due to the finite size of the system seems to be 
lacking. 

In this paper, we discuss finite size effects on the 
dynamics of homogeneous nucleation in a first-order 
temperature-driven transition. In particular, we consider 
the case of cosmological phase transitions in the early uni- 
verse [4], and that of a quark-gluon plasma (QGP) decay 
into hadronic matter in a high-energy heavy ion collision 
[5,6]. The former might provide sensible mechanisms to 
explain the baryon number asymmetry in the universe 
and primordial nucleosynthesis [7,8], whereas the latter 
is expected to be observed [9] at BNL Relativistic Heavy 
Ion Collider (RHIC). The length and time scales involved 
in each of these cases differ by several orders of magni- 



tude. We will argue that, despite statements to the con- 
trary in the literature, the finite size of the early universe 
plays no important role in the process of homogeneous 
nucleation during the quark-hadron and the electroweak 
cosmological phase transitions. For an expanding QGP 
plasma, finite size effects might prove to be relevant if 
the plasma is formed in relativistic heavy ion collisions. 

In the usual description of homogeneous nucleation [1], 
there are two ways in which the finite size of the system 
can affect the formation and evolution of bubbles and, 
consequently, the dynamics of phase conversion. Firstly, 
one has to consider the effects on the nucleation rate and 
the early stage growth of the bubbles. As will be shown 
below, this correction comes about through an intrinsic 
uncertainty in the determination of the supercooling un- 
dergone by the system. For the cases considered here, it 
brings only minor modifications to a description which 
assumes an infinite system. The second and, in general, 
most important finite-size effect is its influence on the 
domain coarsening process, or late-stage growth of the 
bubbles. The relevant length scales here are given by the 
typical size of the system, the radius of the critical bub- 
ble and the correlation length. One has to compare these 
scales, in each case of interest, to measure the relevance 
of possible corrections due to the finite size of the system. 

The arguments presented in this paper are based solely 
on analytic considerations and simple estimates. To ob- 
tain precise results, one should perform numerical sim- 
ulations. Lattice methods have been successfully ap- 
plied to the study of homogeneous nucleation in different 
contexts [3,10,11]. A systematic lattice investigation of 
bubble nucleation and spinodal decomposition, includ- 
ing finite-size effects, in the case of high-energy heavy 
ion collisions is certainly important and will be reported 
in a future publication [12]. 

The paper is organized as follows. In Section II, we 
present a short review of relevant aspects of homogeneous 
nucleation that are important for our analysis. In section 
III, we discuss the effects of the finite size of the system 
in general, and apply the results to the cases mentioned 
above. Section IV contains our conclusions and outlook. 
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II. HOMOGENEOUS NUCLEATION IN A 
FIRST-ORDER PHASE TRANSITION 

In a continuum description of a first-order phase transi- 
tion, it is convenient to make use of a coarse-grained free 
energy, F, which is expected to have the familiar Landau- 
Ginzburg form with temperature-dependent coefficients 
[1]. (In the case of QCD, such a free energy can be ob- 
tained, for instance, from the one-loop effective potential 
of a linear sigma model coupled to quarks [13,14].) For 
temperatures between the critical temperature, T c , and a 
temperature that characterizes the spinodal region, T sp , 
the minimum corresponding to the symmetric phase is 
metastable, and gradually disappears as the system ap- 
proaches T sp , where the symmetric phase becomes un- 
stable. In other words, at T = T sp , the barrier that sep- 
arates the metastable symmetric phase minimum from 
the true vacuum, within this temperature range, van- 
ishes. The driving mechanism for the phase transition 
above the spinodal temperature is the nucleation of large 
localized domains, or bubbles, of the true vacuum phase, 
inside the so-called false vacuum or metastable phase, via 
thermal activation. Below T sp , one has the phenomenon 
of spinodal decomposition, where the transition is driven 
by long- wavelength small fluctuations [1]. An expanding 
system will probe these two regions at a pace given by its 
rate of expansion. The investigation of the time scale for 
thermal nucleation relative to that for the expansion, to 
evaluate the role played by nucleation in the phase con- 
version, is therefore an interesting issue (see, for instance, 
Ref. [14]). However, for a more complete and consistent 
analysis, one also has to incorporate the effects due to 
the finite size of the system under investigation. 

The nucleation rate, namely, the probability per unit 
volume, per unit time, to form a critical bubble, can be 
expressed as T = V e~ F *>/ T , where Ft, is the free energy 
of a critical bubble. The prcfactor V, discussed further 
below, measures statistical and dynamical fluctuations 
about the saddle point of the Euclidean action in func- 
tional space. The critical bubble, with radius R = R c , 
is a field configuration given by a radially symmetric, 
static solution of the Euler-Lagrange field equations cor- 
responding to an exact balance between competing vol- 
ume and surface energy contributions. The critical bub- 
ble is unstable with respect to small changes of its radius; 
for subcritical bubbles, R < R c , the surface energy dom- 
inates, and the bubble shrinks and vanishes, while for 
supercritical bubbles, R > R c , the volume energy domi- 
nates, and the bubble grows driving the decay process. 

It is convenient to write the prefactor V as a product 
of the bubble's growth rate and a factor proportional to 
the ratio of the determinant of the fluctuation operator 
around the bubble configuration relative to that around 
the homogeneous metastable state [15]. For the relativis- 
tic case, in the thin-wall limit, we have [16]: 



16/^3/2 a V R c 

Here, n and £ are respectively the shear viscosity and the 
correlation length in the symmetric phase and Auj is the 
enthalpy density difference between the two phases, a is 
the surface tension of the interface, which is related to 
the critical radius and the difference in pressure between 
the mestastable phase and the true vacuum, Ap, as R c ~ 
2a I Ap. The thin- wall limit is a good approximation in 
the limit of small supercooling, 9 = 1 — T/T c < 1. In 
this limit, the decay rate can be written as 

t, ^ ( 47rer 2 \ _ / 16tt CT 3 \ 
r = ?e X p{-—Rl)=Ve X p{- — , 

(2) 

where I is the latent heat density. For increasing values 
of the supercooling, 9, the argument of the exponential 
decreases and the decay rate increases. Equivalently, the 
critical radius becomes smaller, so that it is easier to 
nucleate supercritical bubbles. The only way in which 
the finite size of the system can affect the decay rate is 
through its influence on the supercooling. (Of course, 
one has to check whether the size of the system allows 
for the presence of at least one critical bubble, including 
the wall that separates the two phases, or not.) 

III. FINITE-SIZE EFFECTS 
A. Coarse-grained free energy 

Before considering the effects of the finite size of the 
system on the nucleation process, let us move one step 
backwards to examine a point that will be relevant for 
the following, namely, the definition of metastable states 
in equilibrium statistical mechanics. Since the complete 
partition function is dominated by configurations that 
minimize the free energy and correspond to equilibrium 
states, it does not allow for the existence of metastable 
branches in the thermodynamic limit. One solution to 
this difficulty consists in constraining the phase space 
and computing restricted partition functions. In order 
to define an optimal cut, one has to resort to physical 
arguments [1]. 

Such a procedure can be implemented as follows [17]. 
The system of size L is divided into cells of linear size 
A cff centered at positions x . The length scale which de- 
fines the coarse-graining of the system, \ cg , should be 
appreciably larger than the underlying lattice spacing a. 
Moreover, within each cell, the relevant order parameter, 
4>, should vary smoothly in space and the equilibration 
time of the system should be much faster than the pro- 
cesses under consideration. This allows us to adopt a 
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continuum order parameter, <p(x,t), denned as an aver- 
age over each cell. 

However, by construction, the coarse-grained free en- 
ergy will depend on the scale X cg and, as mentioned pre- 
viously, one has to appropriately constrain the phase 
space to allow for metastable states. If X cg is small 
enough as compared to the correlation length, £, phase 
separation within a cell cannot occur and one can de- 
fine a coarse-grained free energy, F, such as the one we 
considered in the discussion of Section II. The condi- 
tion \ cg << £ therefore represents the physically moti- 
vated restriction to the partition function that allows for 
a well-defined coarse-grained description of metastable 
states. As a corollary, increasing the value of the coarse- 
graining length relative to the correlation length would 
incorporate more and more fluctuation modes, which 
were excluded by this ultraviolet cut-off. The free en- 
ergy F would then flatten out completely in the limit 
^cg/S, —>■ co, thereby approaching the equilibrium de- 
scription with no metastable branch [18]. For a system 
to be characterized by a well-defined Landau- Ginzburg 
coarse-grained free energy one therefore requires a clean 
separation of scales: a « X cg « £ << L. 

B. Nucleation rate 

The only effect that might be relevant to the nucleation 
rate, T, is the rounding of singularities, since they will 
affect the degree of supercooling. The relevant physical 
quantity can be calculated in a finite-size scaling frame- 
work [19], resulting in the following expression: 

^Fround ^F c 

where L is the typical length scale of the system and 
d is the number of dimensions. The quantity AF roun d 
is a measure of the smoothening of singularities due to 
the finite size of the system [20]. It corresponds to the 
(now non-vanishing) width in temperature of the region 
where the energy density of the system suffers an abrupt 
jump, around T = F c [21]. This phenomenon enters the 
decay rate as a minimum value for 9 which reflects the 
uncertainity in the supercooling due to the finite size of 
the system. This conclusion should be contrasted with 
the arguments presented in Ref. [22], particularly in the 
case of cosmological phase transitions [23] . 

For phase transitions in the early universe, L is given 
by the radius of the universe at a given time (or temper- 
ature). Going back to very early times, t ~ 10~ 4 s and 
T ~ 10 12 K ~ 10 _1 GeV, we can use the elementary 
particle model [24] to obtain an approximate equation of 
state (ideal gas) 

3p«6«^iV(T)T 4 , (4) 



where N(T) = N bosons {T) + {7/8)N fermions (T) is the 
total number of degrees of freedom, e is the energy density 
and p is the pressure. From Einstein's field equations, one 
can derive the relation between time and energy density 
in the early universe, which links the age of the universe 
and its temperature in the following way 

t ~ ( 3 ^ 1/2 ~ -1 (-^_\ 1/2 Mil ^ 

\327rGe) ~ 4tt \ttN(T) J T 2 ' { > 

where Mpi is the Planck mass and G is Newton's grav- 
itational constant. The radius of the universe, as given 
by the particle horizon in a Robertson- Walker spacetime 
[25], dh(t) = t/(l — n), where the scale factor grows as 
a(t) ~ i 2 /3(i+«>) = fn an( j p _ we ^ nag tnc f n ow j n g form: 

L univ (T) « i- (^f)) ^ ■ (6) 

For the assumed equation of state, w = 1/3 and n — 
1/2. If we enter temperature in GeV, the typical length 
scale of the early universe is given, in fcrmi units, by 
L univ {T) « A/T 2 ^fN(T), where A = 1.45 x 10 18 . 

It is clear that, due to the large factor ~ 10 18 (coming 
from Mp{), L un i V will be of importance for supercool- 
ing only for extremely high values of the temperature. 
Since N(T) w 50 for the cases of the QCD and the elec- 
troweak phase transitions, one would need T ~ 10 8 GeV 
to be subject to appreciable finite-size effects on the nu- 
cleation rate. Therefore, since T® CD ~ ©(lO" 1 ) GeV 
and Ff w - O(10 2 ) GeV, these effects are negligible in 
such cases. They might prove to be relevant in the case 
of GUT's transitions, where the critical temperature will 
be much higher. 

For the first order transition of a quark-gluon plasma 
into a hadronic gas in a high-energy heavy ion colli- 
sion, one can easily estimate the role of finite-size ef- 
fects on the nucleation rate. The parameters entering 
Eq. (3) are approximately given by T c « 150 MeV [26], 
I w AB w 4 x (150 MeV) 4 , where B is the bag constant, 
and L w 10 fm. The minimal amount of supercooling 
undergone by the plasma is then negligible, of the or- 
der of 10~ 3 . (For more conservative, smaller values of L, 
the minimum value of 9 is still very small.) The system 
can in principle probe the entire domain in 9, and there 
are no constraints on homogeneous nucleation from this 
side. The question of how fast the system probes the nu- 
cleation region, before reaching the spinodal regime, as 
compared to the nucleation rate was discussed in Ref. 
[14,27,28]. Recently, it has been speculated that the 
RHIC data suggest an explosive hadron production due 
to a rapid variation of the effective potential for QCD 
close to T c . The theoretical reasoning is based on the 
results of the Polyakov Loop Model [29] for the deconfm- 
ing phase transition, which lead to a very fast spinodal 
decomposition regime [30]. 
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C. Late-stage growth 

After the nucleation of a given supercritical bubble, it 
will grow with a certain velocity. The set of all supercrit- 
ical bubbles created integrated over time will eventually 
drive the complete phase conversion in a finite system. 
The scales that determine the importance of finite-size 
effects are the typical linear size of the system, the ra- 
dius of the critical bubble and the correlation length. For 
the reasons discussed in the last subsection, the case of 
cosmological phase transitions is to an excellent approxi- 
mation free from finite-size corrections. The enormous 
numerical value of the Planck mass washes out every 
other scale. Therefore, throughout this subsection we 
will address the quark-gluon-hadron phase transition in 
heavy ion collisions. 

For definiteness, let us assume our system is character- 
ized by a coarse-grained free energy of the form 



F((j>,T) = j d d x 



i(V0) 2 + [/(</>, T) 



(7) 



where U (0, T) is some Landau- Ginzburg potential whose 
coefficients depend on the temperature, and <f>(x, t) is a 
scalar field [31]. For the cases to be considered in this pa- 
per, the order parameter, (f>, is not a conserved quantity, 
and its evolution is given by the time-dependent Landau- 
Ginzburg equation [1] 



dt 



5F 



= 7 [V 2 0-C/'(0,T)] 



(8) 



where 7 is the response coefficient which defines a time 
scale for the system. Equation (8) is a standard reaction- 
diffusion equation, and describes the approach to equi- 
librium [32]. 

If U(4>, T) is such that it allows for the existence of 
bubble solutions (taken to be spherical for simplicity), 
then supercritical (subcritical) bubbles expand (shrink), 
in the thin-wall limit, with the following velocity: 



dR . . 



1 



1 



R(t) 



(9) 



where R c = (d— l)a/AF and AF is the difference in free 
energy between the two phases. Equation (9) is an ex- 
ample of the Allen-Cahn equation [1], which relates the 
velocity of a domain wall to the local curvature. The 
response coefficient, 7, can be related to some character- 
istic collision time as will be done later. 

The description of the late-stage domain coarsening is 
given by the Kolmogorov-Avrami theory [1] , which con- 
tains the following assumptions: (i) bubbles grow with- 
out substantial deformation and are uncorrelated; (ii) 
the nucleation rate is a constant; (iii) the bubble growth 
velocity is constant, v — dR/dt = j(d — l)/R c (limit 
R — > 00). The quantity which is usually computed is the 
volume fraction of the stable phase. 



However, one can measure the importance of finite-size 
effects for the case of heavy-ion collisions by comparing, 
for instance, the asymptotic growth velocity (R » R c ) 
for nucleated hadronic bubbles to the expansion veloc- 
ity of the plasma. In the Bjorken picture, one assumes 
that the central rapidity region exhibits longitudinal ex- 
pansion, so that z(t) = v z t, where v z is the collective 
fluid velocity. Conservation of total entropy leads to 
adiabatic expansion and the following cooling behavior: 
(t/U) = {Ti/Tf. The typical length scale of the expand- 
ing system is then 



L(T) w (v z t c ) 



= L 



(10) 



where L = L(T C ) is the initial linear scale of the system 
for the nucleation process which starts at T < T c . 

The relation between time and temperature provided 
by the cooling law that emerges from the Bjorken picture 
suggests the comparison between the following "veloci- 
ties" : 



Vb 



dR 
dT 



3b£L a \ 
2v z o-T?) 



T 



(11) 



the asymptotic bubble growth "velocity" , and the plasma 
expansion "velocity" 



vl 



dL 
dT 



3£o (Tc 



(12) 



The quantity b is a number of order one to first ap- 
proximation, and comes about in the estimate of the phe- 
nomenological response coefficient -){T) « b/2T (see Ref. 
[33] for details). 

Using the numerical values adopted in the previous 
section and cr/Tj? ~ 0.1, we obtain 



(Ik -I 
v L ~ v z \T 



(13) 



One thus observes that the bubble growth velocity be- 
comes larger than the expansion velocity for a supercool- 
ing of order 9 w v z /2Q < 5%. A simple estimate points 
to a critical radius larger than 1 fm at such values of 
supercooling (see also [14]). Therefore, finite-size effects 
appear to be an important ingredient in the phase conver- 
sion process right from the start in the case of high-energy 
heavy-ion collisions. 



IV. CONCLUSIONS AND OUTLOOK 

We have studied the effects of the finite size of the 
system on the nucleation of bubbles in thermally-driven 
first-order phase transition. The physical situations we 
considered were those of the cosmological transitions in 
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the early universe and the hadronization of a quark-gluon 
plasma after a high-energy heavy ion collision. We have 
shown, on very general grounds, that during the quark- 
hadron and the electroweak cosmological transitions the 
universe is large enough so that one is safe to approximate 
it by an infinite system. Corrections due to its finite 
size are negligible. For the case of heavy ions, although 
the influence on the nucleation rate is also small, effects 
related to domain coarsening and late-stage growth are 
important and should be taken into account in a detailed 
study. 

In order to address the question of hadronization after 
a first-order transition, one should then perform finite- 
size real-time lattice simulations. One would thereby 
avoid the drawbacks implied by analytical approxima- 
tions, such as the thin-wall hypothesis. Moreover, one 
would be able to control the behavior of domains and 
study, for instance, scaling properties [12]. It is also in- 
teresting to study the hydrodynamics of nuclear mat- 
ter at chiral limit as a phcnomcnological description of 
the chiral transition in an expanding quark-gluon plasma 
[34-37]. Results in this direction will be presented else- 
where. 
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